Rectification of Swimming Bacteria and Self Driven Particle Systems by Arrays of 

Asymmetric Barriers 
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We show that the recent experimental observation of the rectification of swimming bacteria in a 
system with an array of asymmetric barriers occurs due to the ballistic component of the bacteria 
trajectories introduced by the bacterial "motor." Each bacterium selects a random direction for 
motion and then moves in this direction for a fixed period of time before randomly changing its 
orientation and moving in a new direction. In the limit where the bacteria undergo only Brownian 
motion on the size scale of the barriers, rectification does not occur. We examine the effects of steric 
interactions between the bacteria and observe a clogging effect upon increasing the bacteria density. 

PACS numbers: 87.17. Jj,05.40.Fb 



OO 



o 



I 

-a 
c 

o 
o 



> 

o 

od 
o 

o 



X 



Recent experiments have shown that swimming bacte- 
ria undergo a rectification effect when placed in a con- 
tainer that has an array of funnel-shaped barriers 
The initial bacteria density p\, is constant throughout 
the container, but over time, bacteria accumulate on one 
side of the container due to interactions with the barri- 
ers. The ratio r = p^ / p^ of the bacteria density p*- 1 -* on 
one side of the container to the density p^ on the other 
side of the container increases over time in the presence 
of a single line of barriers, reaching a saturation value 
of around r = 2.75. By taking advantage of this recti- 
fication effect, it has been experimentally demonstrated 
that various types of bacteria patterns can be induced to 
form, as well as closed circuits such as a bacteria pump 
created from a series of barrier arrays. If dead bacteria or 
genetically modified nonswimming bacteria are placed in 
the system instead of normal swimming bacteria, the rec- 
tification effect vanishes. The dead or nonswimming bac- 
teria still undergo thermal motion and can be regarded 
as Brownian particles. In the experiments, a mixture 
of swimming and nonswimming bacteria was separated 
using the rectification effect of the barrier array. 

The basic questions raised by these experiments arc: 
(i.) what are the crucial properties necessary to give rise 
to rectification of the swimming bacteria and (ii.) what 
prevents the nonswimming bacteria from undergoing rec- 
tification. Possible factors that could be important in- 
clude hydrodynamic interactions with the barrier walls, 
the details of the placement of the driving motor on the 
bacteria, the elongated shape of the bacteria, the distance 
that swimming bacteria move before tumbling and select- 
ing a new swimming direction, and the dimensions of the 
asymmetric barriers. Additionally, other complications 
such as bacterium-bacterium interactions and chemical 
sensing could also play a role in the rectification. In 
Ref. P| , the rectification is attributed to a funnel mecha- 
nism arising from interactions between the bacteria and 
the barriers. 

Motivated by ideas in Ref. |l| , in this work we study a 



simple model for pointlike swimming bacteria which pro- 
duces rectification in the presence of asymmetric barriers. 
Each bacterium moves ballistically under the influence of 
a motor force in a randomly chosen direction for a fixed 
distance lb before selecting a new random direction and 
moving ballistically again. This type of motion is known 
to occur for swimming bacteria which move in a fixed di- 
rection for a period of time before undergoing a tumbling 
locomotion and moving in a new direction [2| . We assume 
that upon interacting with a barrier, the bacterium does 
not reorient but moves along the barrier wall at a ve- 
locity determined by the component of the motor force 
that is parallel to the wall. As a result, the bacteria are 
entrained along the barrier walls as in Ref. [l[ and recti- 
fication of the bacteria density occurs. The rectification 
persists even when an additional Langevin noise term is 
added to the ballistic motion of the bacteria. 

For a system with a single line of barriers, we show that 
the magnitude of the rectification increases with increas- 
ing ballistic flight distance lb, while in the limit where the 
bacteria move like Brownian particles on the length scale 
of the barriers, the rectification vanishes. We also find 
that in this model the bacteria density is highest along 
the container walls, in agreement with experimental ob- 
servations. When steric interactions between the bacteria 
are included, the rectification is reduced due to a clog- 
ging effect. We show that the experimental system is a 
realization of a correlation ratchet Q , in which particles 
with fluctuations that violate detailed balance can give 
rise to a ratchet effect in the presence of an additional 
asymmetry such as an asymmetric substrate. Our model 
should also be useful for other systems of self-propelled 
particles where the particle motion is non-Brownian on a 
certain length scale. Recent examples of such systems in- 
clude self-propelled colloids [1] , artificial swimmers made 
from magnetic beads [B| , autonomously moving catalytic 
nanorods 0], and self-propelled nematics Q. Since the 
amount of rectification is affected by the ballistic flight 
length lb, the asymmetric barriers could be used to sort 
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FIG. 1: Images of the simulation system with 9 = 60°. Black 
dots: bacteria positions; heavy black lines: barrier locations. 

(a) Initial system configuration with uniform bacteria density 
for Nb = 980 and It = 40. The density in the upper chamber 
is denoted p' 1 ' and the density in the lower chamber is p^ 2 \ 

(b) The same system after lOOr simulation steps showing the 
higher bacteria density in the upper chamber, p^ > p^ 2 \ 



different species of swimming or self-propelled particles. 

Our simulation consists of a bounded two-dimensional 
box of an array of funnel shaped barriers as shown in 
Fig. 1. The system has dimensions of L x L, where L = 99 
unless otherwise noted, and contains Nb particles repre- 
senting bacteria at an overall density pb — Nb/L 2 . The 
particles move in an overdamped media where we do not 
take into account any hydrodynamic effects. A single 
bacterium i obeys the following equation of motion: 



(i) 



Here Ri is the position of bacterium i and rj is the phe- 
nomenological damping constant which we set to rj = 1. 
The driving force from the bacterial motor F™(t) has 
a constant magnitude |F m | = 2.0 but is applied in a 
randomly chosen direction which changes after every r 
simulation time steps. In the absence of other forces 
the bacterium would move ballistically over a distance 
lb = TSt\F m \ before reorienting, where St = 0.005 is the 
magnitude of a simulation time step. The resulting ran- 
dom walk motion is illustrated in Fig. 2(a) for r = 1000 
and lb — 10. As r — > 1, the motion approaches a random 
walk even at the smallest size scales, as shown in Fig. 2(b) 
for 4 = 1- The thermal force arising from fluctuations 
in the solvent is given by Ff, which has the properties 
(F7(t)) = and (if \t)Fj \t')) = 2r 1 k B TS ij S{t-t). The 
bacteria motion under the influence of both thermal fluc- 
tuations and the ballistic motor force is illustrated in 
Fig. 2(c) for lb = 10 and F T = 10. For most of this work 
we set F T = 0. The force from the barriers and con- 
tainer walls is given by Ff . Each barrier is modeled by 
two half-parabolic domes of strength fs = 30 and radius 
tb = 0.05 separated by an elongated region of length 
Lb = 5.0 which repels the bacteria perpendicular to the 
trap axis: Ff = £^ [^e( ri )R± + ^efo)^]. 
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FIG. 2: Dots: position of one bacterium in the system. Lines: 
trajectory of the bacterium over a fixed time period, (a) A 
system with only a motor force showing ballistic motion and 
a random reorientation each time the bacterium has moved 
a distance lb = 10. Here F T = 0. (b) A system with small 
lb — 1 and F T — 0. (c) A system with both ballistic motion 
and thermal motion, where lb = 10 and F T — 10. 



Here the total number of barriers (including confining 



walls) is Ni 
Rt 



IR, 



= 28, n 
-RfiiBPffl 



= ri 
R ife 



R 



± 



R ik = 



■tk> r2 = TB 

(Ri - Rf ± £bP|| )/ R ik , 



(Ri - Rf ) • pi|, Ri = ((Ri - Rf ) • P k ± )/Ri, 
Ri (Rf ) is the position of bacterium i (barrier fc), and 
p| (p*l) i s a unit vector parallel (perpendicular) to the 
axis of barrier k. Each funnel is modeled as two barri- 
ers meeting at a common endpoint and placed at angles 
9 and ir — 9 with the a;-axis. Our system contains 12 
funnels in a one-dimensional array with a lattice con- 
stant of Is = 8.25. The steric force from the bacterium- 
bacterium interaction, Ff , is modeled as a stiff repul- 
sive spring of range r s = 0.35 and strength f s = 150, 
with Ff = V^,^Bir : .iR, ; . Here r 3 = r s - R %] \ 

Rij = | R ? ; — Rj | , and Rij = (Ri — Rj ) / R%j ■ For 
most of this work we consider the limit where bacterium- 
bacterium interactions are irrelevant and set f s = 0. 

We base our system geometry on the work in Ref. [l[. 
In Fig. 1(a) we illustrate the system geometry, showing 
the array of funnel shapes, the boundary walls, and the 
initial locations of the bacteria. The bacteria density in 
the upper (lower) chamber is /r (p( 2 )). The angle 9 of 
each funnel barrier is set to 9 = 60°. We first investigate 
the simplest set of parameters that produce rectification 
by fixing F T — 0, f s — 0, and = 40. The initial density 
ratio r — p^ / p^ = 1, as shown in Fig. 1(a). In Fig. 1(b) 
we illustrate the same system after 100r time steps when 
r = 2.3, showing a density increase in the upper half 
of the sample. A buildup of bacteria occurs along the 
barrier walls and the boundaries, similar to that observed 
in the experiments of Ref. [l| . 

In Fig. 3(a) we plot the time dependent behavior of 
p^ / Pb and p^/pb for the same system shown in Fig. 1 
with lb = 180. We also show the same quantities for a 
system with a much smaller lb = 0.01 in which the bac- 
teria motion is equivalent to a random walk on the size 
scale of the barriers. For the small lb there is no rectifica- 
tion of the bacteria density as would be expected, while 
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FIG. 3: (a) The time dependence of the bacteria density 
P I Pb (upper line) and p' 2 ' / p b (lower line), normalized by 
the overall bacteria density pb, for the same system in Fig. 1 
with lb — 180, showing the density buildup over time in the 
top half of the sample. Dotted lines: p^ / pb an d p^ / pb 
for the same system but with lb = 0.01 so that the bacteria 
undergo Brownian motion on the size scale of the barrier and 
no rectification occurs, (b) r = p'^/p' 2 ' versus h for a sample 
with = 60° measured after 10 3 r simulation steps. 

for the large 4 the system rectifies. The time dependence 
of r = p^'/p^ 2 ' at short times can be fit to an exponen- 
tial form, similar to the experiments of Ref. We note 
that asymptotically, r 4.3, while the asymptotic value 
of r depends on 4- We have repeated the simulations for 
different values of 4, and in Fig. 3(b) we plot r versus 
4 measured after 10 3 r simulation time steps. For low 
lb < 200, the rectification increases monotonically with 
lb, with r reaching a maximum of value of r = 4.8 at 
lb = 200. We note that r depends in a similar way upon 
the barrier length Lb for fixed lb, and the rectification 
vanishes for large values of Lb when there is no longer a 
gap between adjacent funnel barriers. 

These results suggest that hydrodynamic effects or 
chemotaxis are not required to produce the rectification 
seen in Ref. but that instead the rectification arises 
from the fact that the bacteria move in finite steplike 
motions. In Fig. 4(a) we illustrate the specific mechanics 
that give rise to the funneling behavior of the bacteria at 
the barrier. As described in Ref. a bacterium which 
collides with a barrier cannot continue to move in the di- 
rection of the motor force; instead, the bacterium slides 
along the barrier at a velocity determined by the compo- 
nent of the motor force parallel to the barrier. If lb is large 
enough, there is a high probability that the bacterium 
will slide all the way across the barrier to the other half of 
the container before the motor force changes its orienta- 
tion. If lb is very small, the bacterium cannot move very 
far along the barrier before it randomly moves off when 
the motor force changes, as illustrated by the light line 
in Fig. 4(a). The rectification is produced by the combi- 
nation of the asymmetric barrier and the motion of the 
bacteria that allows a breaking of detailed balance when 
the bacteria interact with the barriers. Another conse- 
quence of this model is that the bacteria should tend to 
accumulate on the container barrier walls, as observed in 
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FIG. 4: (a) Schematic trajectories for: (dark line) a bac- 
terium moving under motor forces only and being entrained 
by the wall; (light line) a bacterium experiencing only Brow- 
nian forces. The barrier spacing Is, opening size l , funnel 
width l w , and angle cj> between the bacterium trajectory and 
the funnel barrier wall are also indicated, (b) r = p^/p^ 
vs for a system with lb = 20, F T — 0, f s — 0, and L — 99. 
Inset: r vs for the same system with fixed l = 3.25 and 
varied L. (c) r vs F T for a system with lb = 20, 8 = 60°, 
L = 99, and f s = 0. Inset: r vs Lb for a system with l\, — 20, 
F T = 0, f s — 0, fixed l = 3.25, and varied L. (d) r vs l for 
a system with lb = 20, 9 — 60°, Lb = 5.0, f s ~ 0, and varied 
L. Inset: r vs overall bacteria density pb for a system with 
steric interactions, f s = 150, 6 — 60°, f T = 0, L = 99, and 
h = 20. 



the experiment [lj. We note that the real bacteria may 
tend to be hydrodynamically trapped by the wall. This 
effect would tend to enhance the rectification effect since 
the bacteria could move for distances longer than lb along 
the barrier, generating rectification even for moderately 
small 4- 

To compare our simulation results quantitatively to ex- 
periment, we take a simulation distance unit of 1 micron 
and time step unit of 5 x 10 -4 sec. This gives a bacteria 
velocity of 20 microns/sec, close to the estimated speed 
v « 23 microns/ sec for E. coli [1, Q . A tumbling rate of 
roughly 1 Hz 0] corresponds to 4 = 20 microns; Ref. [l[ 
gives the estimate lb ~ 50 microns. Fig. 3(a) shows that 
our system reaches a stationary state within 1 x 10 5 time 
steps, which corresponds to 50 minutes. This compares 
well with the saturation time scale in the experiment of 
20 to 80 minutes 

We next consider changing the geometry of the barrier 
funnels. In Fig. 4(b) we plot r vs 0, the funnel angle, after 
10 3 simulation steps for a system with 4 = 20 and fixed 
L = 99. For 8 < 41° adjacent funnel barriers are in con- 
tact so there is no transport. The maximum of r = 4.6 
occurs at — 41°, and r decreases with increasing 9 for 
6 > 41°, reaching r = 1 at 8 = 90°. The experiments also 



4 



showed no rectification for fiat barriers The inset of 
Fig. 4(b) shows the dependence of r on 9 when the space 
between funnel barriers l a is held fixed, indicating that 
the maximum rectification occurs for 9 = 60°. The inset 
of Fig. 4(c) illustrates that r initially increases with bar- 
rier length Lb before saturating, while Fig. 4(d) shows 
that r decreases as l increases. 

A simple heuristic expression gives a lower bound for 
the value of r in the case when the spacing Is between 
adjacent funnel barriers is much smaller than the ballistic 
flight length, 1$ <C h- Assuming that the bacteria are 
entrained by any barrier wall which they strike, detailed 
balance implies that in the steady state at long times, 

r = P [1) /p {2) = Prob.{(j> > n/2)(l w /l ) + 1. (2) 

Here, the tip-to-tip width of a single funnel barrier is l w = 
Is — l , as indicated in Fig. 4(a). In Eq.flSJ, Prob.(4> > 
7r/2) is the probability that whenever the bacteria hit 
the funnel barrier walls, they do so at a large enough 
angle <f> such that they slide along the wall in an upward 
direction to emerge in the top portion of the sample, as 
shown in Fig. 4(a). Equation ([2|) leads to a value of r 
which is strictly bounded from below by one, reflecting 
the rectification effect. By measuring Prob.(4> > 7r/2) we 
have verified that Eq.((2]) is indeed generally satisfied. 

In the experimental system, the bacteria experience 
some additional Brownian fluctuations due to the sol- 
vent. We model this effect by setting F T > 0, pro- 
ducing the type of motion shown in Fig. 1(c). If only 
the thermal noise term is applied and the motor force 
F M = 0, we find no rectification. In Fig. 4(c) we show 
r vs F T for a system with lb — 20 and barrier radius 
T6 = 0.5. As F T increases, the rectification is slowly 
reduced, while for large enough F T the rectification is 
washed out and r = 1. The destruction of the rectifi- 
cation occurs because the thermal term causes some of 
the bacteria to move away from the barriers rather than 
channeling along them, so the progress along the walls is 
reduced from the mechanism shown in Fig. 4(a). 

The bacteria in the experimental system have a steric 
interaction with each other which prevents overlap. To 
mimic this effect we add a steric force interaction by set- 
ting f s > 0. In Fig. 4(d) we plot r vs bacteria density 
Pb = Nb/L 2 for the same system in Fig. 2 but with a 
finite steric interaction strength of f s = 150. As the den- 
sity pb increases, the effective rectification is reduced due 
to a clogging effect since only a small number of bacteria 
can fit along the barrier wall. The effect of the clogging 
can be reduced by the addition of thermal noise since this 
has a tendency to prevent bacteria from clustering. 

We have also performed simulations for mixtures of in- 
teracting bacteria in which one species moves under both 
a motor force and thermal noise while the other species 
moves only due to thermal noise. This mimics a sys- 
tem with a combination of swimming and nonswimming 
bacteria. If the density of the nonswimming bacteria is 



sufficiently high, the rectification efficiency is reduced. It 
is also possible for an entrainment effect to occur in which 
the swimming bacteria can induce a weak rectification of 
the nonswimming bacteria. 

The rectified motion in this system has similarities to 
Brownian ratchets where particles undergoing Brownian 
motion in the presence of an asymmetric potential can 
exhibit a net drift in the presence of an additional ac 
drive or flashing potential substrate [X0| . Our system is 
not a Brownian ratchet in this sense; instead, it can be 
considered to be a realization of a correlation ratchet. In 
correlation ratchets, an overdamped particle can exhibit 
dc drift on an asymmetric substrate in the absence of 
an ac flashing or rocking provided that the fluctuations 
of the particle motion have certain properties that break 
detailed balance [3|. In bacteria and other self-propelled 
particle systems, it is the motor force that causes the 
fluctuations to break detailed balance. 

To summarize, we show that a rectification phe- 
nomenon similar to that observed in recent experiments 
[H can be achieved in a simple model for overdamped 
swimming bacteria where the bacteria move in straight 
lines of length lb before randomly changing direction 
while also interacting with funnel shaped barriers. In 
this model, no hydrodynamic interactions are needed to 
produce the rectification. In the limit of bacteria which 
undergo only Brownian motion on the size scale of the 
barriers, no rectification occurs, in agreement with ex- 
periment [H- The rectification results when If, is large 
enough that bacteria which strike a funnel barrier move 
along the entire length of the barrier before the bacterial 
motor force changes direction. We have also examined 
the effects of steric interactions on the rectification and 
find a clogging phenomenon at the barrier funnels which 
reduces the effectiveness of the rectification. Addition- 
ally, if a sufficiently strong Langevin noise term is added 
to the model, the rectification is reduced. This system 
can be regarded as a two-dimensional realization of a cor- 
relation ratchet in which fluctuations that break detailed 
balance can give rise to directed motion in the presence of 
an asymmetric substrate. Our results should also be rel- 
evant to other two-dimensional systems of self-propelled 
particles that undergo a similar type of motion. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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